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Abstract 

We first establish the presence of a diffractive front in the fundamen- 
tal solution of the wave operator with a diract delta intial condition in two 
dimensional euclidean space caused by the potentials perturbation on the 
spherical laplacian. This motivates a result which restricts the propagation 
of singularities for the wave operator with a more general potential to pre- 
cisely these diffractive fronts higher dimensional euclidean spaces. This is 
proven using microlocal energy estimates. 

1 Introduction 

This paper deals with the diffractive phenomenon exhibited by the wave equation 
when singularities in its distributional solution encounter an inverse square poten- 
tial. Potentials of this order are of broad interest since they form the border line 
case for the existence of global in time estimates for the wave and Schrodinger 
equations (see HG-VI HR-SIO ; recent work has been done by Burq, Planchon, 
Stalker, and Tahvildar-Zadeh to establish Strichartz estimates for these equations 
with such a potential in IB-P -S-TZl . Inverse square potentials also appear in phys- 
ical equations, including a recast version of the Dirac equation with Coulomb po- 
tential (see [CaJ) and the linearized perturbations of the Schwarzschild solution 
among others for the Einstein equations (see HR-WH 0). Furthermore, the inter- 
action of singularities in solutions of PDE with singular spaces have long been 
known to cause diffraction; we will survey the work done in this direction before 
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stating some diffractive results for inverse square potentials. The mathematical in- 
terest of this potential in all of these cases largely comes from the —2 radial order 
of homogeneity common to both the inverse square potential and the Laplacian. 

In Euclidean spacetime, R'™ -1 x R = {(x,t)}, the simplest version of our 
scenario is the wave equation with an inverse square potential ^ (where r is radial 
distance in x from x — 0) which is described by: 

Du(x, t) + ^u(x, t) = (a > 0, a constant) 

u t (x, 0) = 8(x-x ) (1) 
u(x,0) = 

The "box" symbol: □ is the d'Alembert operator defined as: 

D : ~ at 2 2s dxf { ) 

i=l 1 

where the second term is the familiar Laplacian A := — ~& m me n ~ 1 

spatial variables. 5(x — xq) is a Dirac delta function giving an initial condition with 
its singular peak centered at x . Our interest is in finding the set of singularities of 
the distributional (weak) solution u in space-time. This set is denoted singsupp(w) 
and defined as the places where u is not C°°. 

To illustrate, consider figure 1 below which qualitatively shows the behavior of 
the solution u in the spacetime R 2 x R. In the forward time direction, the behavior 
of the initial data 5(x — xq) is like an "explosion" at t = which sends out a 
primary spherical front of singularities at unit speed that sweeps out a cone (with 
profile slope 1) in space-time, accounting for the outer cone which is between 
regions I and II in the figure; this behaviour is the same as what would occur for 
the free wave equation, Du = 0, under the same initial conditions. The inner cone 
between regions II and III is a diffractive front in the solution brought about by 
the initial front striking the potential at the origin; this is a qualitatively different 
phenomenon from the singularities in the free case. 

A similar diffractive front is the subject of several other results. The first 
rigorous treatment of diffraction for the wave equation was made by Sommerfeld 
in [[Soil for the case of a singular front encountering a slit in a plane, causing 
diffraction into the shadow region of the slit. Friedlander extended this to cone 
obstacles in [FrJ followed by Cheeger and Taylor who explicitly dealt with the 
case of general product cones in HC-TII . The potential can be also thought of as 
a singular feature of our space at the origin, and by employing the separation of 
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Figure 1: main front (with embedded bicharacteristic) and diffractive front. 

variable methods used by the authors above, we can deduce the following theorem 
about the explicit solution: 

Theorem 1.1. Existence of diffraction: There is an explicit diffractive front cor- 
responding to the inner cone in figure 1 in u, the explicit solution o/di]) in IR 2 . 

In more general geometric contexts, we can't use separation of variables to 
achieve a similar explicit solution; energy methods are preferred. Work by Lebeau 
gives a diffractive result in the analytic setting for a broad class of manifolds in 
(|L|. Melrose and Wunsch have established results for conic manifolds in HM-WD 
in the smooth setting, and along with Vasy who studied manifolds with corners 
in the smooth setting [V|, they have recently published results in the case of edge 
manifolds HM- V-WD . Using (microlocal) energy techniques employed by these 
authors, we are led to a result about the propagation of smoothness around points 
of the form (r = 0, t) where the potential is concentrated: 

Theorem 1.2. Propagation of smoothness: Let ube a solution of 

| Uu{x,t) + ^P-U(x,t) = /GC <oO ; || / || oo<00)/ >_(n_2) 

(3) 

If there are no singularities of u going into the point (0,t) originating from 
this point's backward cone of influence, then there are no singularities coming out 
of its forward diffractive cone as seen in figure^ 
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Figure 2: incoming smoothness implies outgoing smoothness. 

Note that there are no other possible incoming or outgoing singularities at this 
point except for those in the two cones depicted since they must all travel with 
unit speed due to a standard result from microlocal analysis. Also note that the in 
Theorem ll.li the single singularity traveling along the bicharacteristic toward the 
origin is enough to cause the shower of singularities that make up the diffractive 
front. Since energy estimates tell us about the propagation of regularity, we can't 
hope for a better result in this context than the propagation of smoothness along 
the light cone out of the singularity given smoothness along the entire incoming 
cone. Results of this type can be thought of as a microlocal version of the local 
smoothing estimates for the wave equation seen in HB-P-S-TZI that take a further 
phase variable into account. 

Because of the arbitrary bounded smooth function in the potential term, we 
can't separate variables to find an explicit solution which directly exhibits diffrac- 
tive behavior. Energy methods become our primary tool, and the resulting theo- 
rems deal with the propagation of smoothness, which in turn restricts the behavior 
of singularities. The result above restricts the outcome of interactions between 
singularities from the solution and the potential's singular point to only the emis- 
sion of a diffractive front immediately following an incoming singularity; we will 
not see the spontaneous emission of singularities from the potential, nor will sin- 
gularties be trapped and released at a later time. 
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2 First Diffractive Theorem 



We first prove that diffraction exists: 



Theorem 2.1. Existence of diffraction: There is an explicit diffractive front cor- 
responding to the inner cone in figure 1 in u, the explicit solution of 



in 



Du(x,t) + $u(x,t) = 

u t {x,0) = S(x-x ) 

u(x,0) = 

. space-time. 



(a > 0, a constant) 



(4) 



Proof. We follow the general method used in HC-TH . Recall that the solution op- 
erator to this perturbed wave equation is (— A)~ 2 sin t(— A) 2 where A = A — -%, 
a E M + , so the fundamental solution is (— A)~5 sint(— A)^5. We would like an 
expression for the fundamental solution in terms of polar coordinates to exploit 
the radial symmetry of the solution, and isolate an expression for the diffractive 
front in terms of the spectrum of A. 

We consider i 2 = R x S 1 as a cone with base S 1 , i.e. in polar coordinates 
(r, 9). (Our 9 E [0, 2n] as usual.) Here, the Laplacian in polar coordinates has the 
form: 

d 2 Id Id 2 

+ -— + -— (5) 



r dr r 2 89 2 
-£r + a) 2 and u j 



(jjj + a) 2 where {/i,} 



dr 2 

Define the operator v on S 1 by v 
is the spectrum for with eigenfunctions <p{9). The base S 1 = W/2nZ, and we 
get v n = {n 2 + a) 2 for n G Z and </?(0) = e me . Separating variables, we consider 
the action of A on g(r, 9) = gi(r)ipj(9) namely: 



A(/(r,i 



E 



9r 2 



r dr 



9j(r) <Pj(0) 



(6) 



The form of the radial part strongly suggests Bessel's equation: 

Ql 1 A ,,2" 



1 V 

dz 2 z dz z 2 
which are solved by the Bessel Functions: 

(z/2Y 



Ju(z) 



-JJz) 



J„ 



r(i/2)r> + 1/2) 



(1 - t 2 y- 1,2 e izt dt 



(7) 



(8) 
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This suggests making use of convolution with the Bessel Functions in a spec- 
tral transformation known as the Hankel Transform: 

H v (g)(X) = [ g(r)J u {Xr)rdr (9) 

for which we know the following proposition: 

Proposition 2.2. For v > H v extends uniquely from to 

H u : L 2 (M+, rdr) -> L 2 (M + , XdX) (10) 
and for each g G L 2 (M + ,rdr): 

H u oH u g = g (11) 

(See IfCTl for a proof.) 

Letting + (£)] , the radial part of the operator we have 

W(MXr)) = - 7 2 ^(Ar) that: 

/■oo 

H u {L u g) = / L u (J u (Xr))grdr (12) 
Jo 



= -A 2 / J u (Xr)grdr (13) 

= -A 2 ^(^) (14) 

and this isometry carries the radial portion of A into multiplication by —A 2 . 

Combining this with the Fourier Transform on the spherical variables gives as 
the Schwarz kernel for the solution: 

POO 

f{-h){r l ,r 2 ,e l ,e 2 ) = Y j e in ^-^ / f (X 2 )J Un {Xr 1 )J Un {Xr 2 )XdX (15) 

we can write the Schwartz Kernel for our operator v for a single Fourier mode v n 
as: 

/■oo 

f(-A)(n,r 2 ,v n )= /(A 2 )J„„(Ar 1 )J,„(Ar 2 )AdA (16) 
Jo 

Apply this to our actual fundamental solution, ^/Xsm{t^/X) = lim e ^ im { ^ lt)x 
and then use the Lipschitz Hankel Integral (see HC-T1I ) to obtain the following 
equality: 

v 1 / r 2 _i_ r 2 _i_ ^2 



e- tA J 1 ^(riA)J I ^(r a A)dA= -(n^)- 1/J g^ ( 1 ^ ) (17) 



where <5^_i is a Legendre function of second kind. 

We take the real solution as the limit of the meromorphic function Q v _x(Z) 
as the Z variable approaches the real line. We will see that Z > 1 corresponds 
to the interior of region I, — 1 < Z < 1 the interior of region II, and Z < — 1 
the interior region III. Approaching the poles at Z = —1,1 from different sides 
will correspond to the boundary of the three regions defined by the cones in our 
diagram and we will show there is a non-zero difference along the entirety of the 
boundary between regions II and III. 

To see this, IlLell gives that Q v _\ has the following integral representation: 



The Legendre function is analytic when Z ^ C\(— oo, 1]; there is a branch cut 
on the real line and analytic continuation allows us to extend both sides of the 
equation to the complement of the cut by shifting our contour of integration on 
the right hand side. As Z — > 1 or Z — > — 1, the integral becomes singular, and the 
directions in which Z approaches these singularities will account for the behavior 
of the solution on different regions of the lightcone. 
Substituting Z = cosher/) yields: 



and note that the cut on C will change under the cosh - transformation in the 
imaginary line. We have coshr/ 6 [—1, 1] <^ cosir] E [—1, 1] <^> ir] E [0, — n] <^ 
i] E [0, in], where correspondences are Z = — 1 to i] = iri, and Z = 1 to i] = 0. 
So the cut becomes [0, iri} U [iri, oo + ni]. The upper half plane is transformed into 
the rectangle enclose by the cut and real line. 

We then have three regions to consider on R: I: cosh(?7) E (l,oo);II: cosh(?7) E 
(— 1, 1); III: cosh(?7) E (— oo, —1) for Q u _i (cosher?)). The contours taken are the 
following: 

In I, we have that the integral is on a purely real contour, and since we are 
concerned with only the imaginary part of Q u _ i (cosh(Ty)), we have K(r 1 ,r 2j v) = 
0. 

In II, we take the contour (see above II) on the cut. As above, the integral 
disappears when the contour is on the real line, leaving only the portion from 
the imaginary segment in [0, iir). We translate this into a real integral on a seg- 
ment between [0, ir). Taking the imaginary part of the result noting the identity 
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(18) 




(19) 
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Figure 3: Equivalent contours for case III, we choose the former for its similarity 
to the contour from case II as we approach ni along [iri, iri + oo]. 



cosh(z) = cos(iz), and that an additional i comes from swapping the order of the 
denominator gives yields: 



,-1/2 



2 , 2 t 2 
-1 / r l+ r 2~' 



COS Vr,S 



it 2 — r\ — r 2 + 1r\r 2 cos s) 1 / 2 



ds 



K(r 1 ,r 2 ,u n ) = -(nr 2 

71 

(20) 

In III, we have several options for contours that go to oo without crossing the cut 
that are equivalent because of analyticity within the rectangle which represents 
the upper half plane under the cosh -1 transform. We choose our contour as an 
extension of the contour in II to facilitate calculating the difference between the 
two (see figure). We add to the portion on [0, in] an additional segment from 7] to 
in and hence an additional term which we again shift into a real integral where 



K(r 1 ,r 2 ,u n ) 



(nr 2 ; 



-1/2 



n 

--{r x r 2 ) 
n 



cos u„s 



(t 2 



r\ + 2r\r 2 cos s) 1 / 2 



ds 



-1/2 



Sill 7TK 



n J 



(2 cosh p - 2 cosh s) 1 !- 



■ ds 



(21) 
(22) 



1 /^.2_i_y,2 £2 \ 

In the integral for case I, we have cosh 7] > 1 so cosh ( 1 2r ^ ) > which 

implies < t < \r\ — r 2 \. We are outside the primary lightcone in case I; the fact 
that our solution is here is consistent with finite speed of propagation. 
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Figure 4: Main front and diffractive front 



In the integral for II, we have " — 1 < cosh 77 < 1" which translates on our 

_ 1 / r 2, r 2_ i 2\ 

contour into < cos ( 1 2r ^ r2 ) < n an d hence \r 1 — r 2 \ < t < r± + r 2 . We 
are inside the primary lightcone but outside of the inner lightcone in case II. 
Finally for case III, "cosh < —1" which really means we approach the cut 

I / £2 ^2 ^2 \ 

from above. We have cosh ( 1 2 1 > and t > r-y + r 2 . This means we are 
inside the inner lightcone. 

Letting r] approach the boundary of case II and III from opposite directions, 
that is ri \ t — r 2 and n f t — r 2 respectively, we see that the solutions in case 
II and III differ by: 

i rP p-su n 

- ~(r^M^) I (2cosh/3 _ 2coshs)1/ ^ <*) 

If we take a second order approximation for cosh near 0, the limit as we ap- 
proach the boundary cone of the the integral above is: 

f 3 e ~ SVn , f 3 ds 

llm / 77; TT> 7, \ — 7777 « s = llm / 7775 77777 ds ( 24 ) 



(25) 
(26) 



p^oJ (2cosh/5-2coshs) 1 /2 ^ J Q (/p _ s 2)i/2 

. s x & 

lim sin 1 



71 

2 
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and so the solution exhibits a difference of: 

-^(r 1 r 2 )- 1 / 2 sin(7r i /„) (27) 

on the n-th mode of the Fourier decompsition of (p(9). 

Recall the spectrum of our base, v n = [n 2 + a) 2 , n E Z. This means if we 
have a 7^ 2nm + m 2 , m E Z, then shi7rz/ n 7^ 0. By uniqueness of the spectrum, 
our solution along the inner cone is not 0. 

We need to eliminate the further possibility of the solution being a collection 
of singularities at discrete points instead of a legitimate front spread across the 
diffractive cone. The Schwarz representation theorem tells us that this undesirable 
case can only come from having a sum of Dirac delta functions and its derivatives 
Y2 k c k S^ k \ which has the spectrum 2~2k>o Cknk - Our solution has the spectrum 
C sin(7rv / n 2 + a) for C some constant; the expression grows asymptotically of 
order 1 so it can only possibly match a sum Dirac delta functions without deriva- 
tives. However, this has the spectrum J2 k>0 c k which is invariant in n and cannot 
match all the values of C sin^V^ 2 + a) over n E N for any values of a 7^ since 
it becomes a-periodic because of the irrational values that \Jn 2 + a takes. Hence 
the fundamental solution on the diffractive cone is a legitimate, non-zero front by 
the uniqueness of spectrum for distributions. □ 

3 Second Diffractive Theorem 

By the first theorem, diffraction indeed occurs when a singularity strikes the ori- 
gin; we would now like to show theorems in more general contexts which restrict 
where singularities may propagate. 

3.1 First Formulation 

Recall that by the celebrated Duistermaat-Hormander theorem, we have away 
from r = 0: 

Theorem 3.1. (Duistermaat-Hormander). Let □ + = 0. WF(u) is a union of 
maximally extended null-bicharacteristics of the vector field H a ( a+ j^) = 2rd t — 

This says the singularities of □ + -% propagate along unit speed geodesies that 
fan out from the singular concentration at x of our initial condition. (When r — > 
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however, we can't propagate past the large singularity.) This theorem is secretly a 
statement about the restriction of the propagation of smoothness in wave operator 
solutions along null-bicharacteri sties of the Hamilton vector field. In the potential 
free case, it in turn restricts the propagation of singularities. 

This means that the single geodesic striking the origin in our first diffractive 
theorem is the sole cause of the entire shower of singularities in the diffractive 
front. Thus, in formulating an analogous propagation of smoothness theorem as 
Duistermaat-Hormander for the wave operator with potentia, we cannot localize 
as finely as we could before since there is no hope of restricting the incoming 
smoothness or the out going smoothness in anything less than an entire cone. 

On the other hand, energy methods allow us to gain some generality in the 
potential: 

Theorem 3.2. Propagation of smoothness: On M n /or n > 3 and \(n) = ^rp, let 
u EV be a solution of 

fir 6) 

Du(x, t) + J -±^u(x, t) = (28) 

where f(r, 6) E C°°{R n / {0}), / > -A(n) 2 locally, and ||/|U,ioc < oo. 

If there are no singularities of u going into the point (0,t) originating from 
this point's backward cone of influence, then there are no singularities coming out 
of its forward diffractive cone as seen in figured 

In the potential, / does not necessarily have a form fi(r)f 2 (9) that separates 
r from 9, so we can't separate variables as before to find an explicit solution and 
resort to energy methods. In particular, we will use microlocal energy estimates 
since boundedness and smoothness of / allow commutator arguments to be made 
in the style of Duistermaat-Hormander and borrow many of the finer techniques 
used in HM-V-WH . We are not getting as finely localized a result as we could on 
our usual Euclidean space, but this result will imply that we don't have sponta- 
neous emission of singularities from the origin due to the potential or trapping of 
incoming singularities with a delayed release. This is an analogous microlocal 
version of the local smoothing estimates that appear in HB-P-S-TZH . 

Our solutions all belong to a domain which on the spatial variables we pick 
to be the Friedrichs extension for the quadratic form associated with our operator 
A+ " 



f_M) . 



Q{u) := / |Vn(a;)| 2 + f -^l\ u (x)\ 2 dx (29) 
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We will have a more thorough discussion of domains after we define the b-Calculus 
and its associated b-Wavefront set as well as b-Sobolev spaces. Ultimately, we 
will want to define all of these objects with respect to our domain, which we will 
prove to be equivalent to b-Sobolev Spaces using Hardy's Inequality. 

We have suppressed the microlocal details of our thereom in its first formula- 
tion above. Indeed, our statement refers to incoming and outgoing singularities, 
which is best formulated in terms of the wavefront set for our solution. This is 
analogous to the following formulation of Duistermaat-Hormander where c(t) is 
a parametrization of our null-bicharacteristic: 

Theorem I3.1P Under the hypotheses o f\3.1\ and WLOG letting our point of 
interest be t = 0, jfWF(u) n {c(t) \ t < — e} = for some e > then 
c(0) i WF(«). 

Due to the prescence of a radial point in the vector field, it will turn out that the 
best propagating variable to use will be the momemtum variable £ corresponding 
to r. However, along constant speed geodesies in Euclidean space, this variable 
makes an abrupt jump as we cross from negative to positive for any geodesic 
striking the origin. It is preferable to use the rescaled rd r to generate our vector 
field and correspond this to the momentum £. This yields a gentle propagation 
through the origin with the momentum vanishing to as it approaches the origin 
and then gradually increasing in speed as it leaves. 

The presence of rd r motivates the use of the b-Pseudodifferential Calculus 
when formulating our commutant and a blowup at the origin in order to accomo- 
date these operators. These can be thought of as the pseudo-differential calculus 
generalized from the vector fields generated by rd r , dg,d t . 

3.2 Blowup and b-Pseudodifferential Calculus 

Here, the result of our blow-up on M n is denoted [M n , 0] = R + x S n-1 ; we are 
continuing polar coordinates from W 1 — {0} in a one-to-one way to the origin by 
joining a sphere. Letting r be our boundary defining function and 6>j represent 
spherical variables, the b- vector fields are generated by rd r , de v d t ; these vectors 
are tangent to the boundary sphere of our blown up space. In particular, note that 
a flow in the rd r direction will slow down as we approach the origin (see figure). 
The bi-filtered *-algebra of pseudodifferential operators associated to these vector 
fields is known as the b-calculus. 

The b-Calculus on a manifold with boundary M, whose formalization is due 
to Melrose (see [MJ) is a microlocal generalization of vector fields that are tangent 
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Figure 5: Blowup with circle at r = 0, rescaled rd r vector field, and propagating 
variable. 

to our blown up boundary, spanned in C°° by rd r , d t and d Zi which correspond to 
the symbols £, r, and Q. The dual to such tangent vectors is b-cotangent bundle 
b T*M. The dual of rd r is — , 8q. is d6i, and d t is dt. We denote the coordinates 
with respect to this basis by £, (, and r thus making our canonical one-form 

+ (idP + rdt. The co-sphere bundle is b S*M = ( b T*M - {0}) /R+ and it 
is here where our characteristic set lives and where our propagation takes place. 

The b-Calculus of pseudo differential operators is denoted \I/™''(M). It is a bi- 
filtered (m for differential order, I for weight order) *-algebra of operators acting 
on C°° and C~°°, smooth functions that vanish to infinite order at the boundary 
and their corresponding dual respectively. They satisfy: 

(I) DifC(M) C y™>°, 

(II) for x a boundary defining function of M, we have x l e *&®' l (M) and 
x l ^'°(M) = ^™\M), 

(III) there is a principal symbol map is a bi-filtered *-algebra homomorphism 
which maps to weighted polyhomogeneous symbols 

V< : K'\M) - x l SZ m ( b T*M\0) (30) 

(IV) the principal symbol sequence is exact 

-> ^~ U (M) ^ *^(Af) -> x l SZ m ( b T*M\0) (31) 



13 



(V) for A e V™' l (M), B e #™''''(M) with principal symbols a and b of respec- 

tive orders, [A, B] e #™+™'-W(m) obeys 

b a m+m ,_ hl+l ,([A,B]) = -{a,b} (32) 

(VI) all A E ^ ' on C°°(M) extend by continuity to a bounded L 2 operator. 

Notice that the original weight of r~ l remains in the residual terms quotiented 
away in the principal symbol map. For our commutator argument which iterates 
on successively on the level of the principal symbol, this will prove to be too 
crude. We will eventually prove lemmas that allow us to explicitly place weights 
on specific differential operators which we factor out at each successive step. 

We say A e is elliptic at p e b S*M if its principal symbol has an inverse 
near p in r-- l C( b S*M) := r- l S h ^ l n ( b T*M)/S h ^ 1 i b T*M). We call these points 
Ell b (A). There is a corresponding concept of microsupport as well, given by the 
essential support of an operator's (left) symbol, denoted WF' b (A). This obeys 

(I) WF\{AB) C WF' b (A)n WF\(B) 

(II) For A e *™''(M) and p e Ell h (A), there is a parametrix Q e ^~ m ^ l {M) 
such that 

p i WP b {QA - I) U WP^AQ - I) (33) 

(III) For A e *&™ ,l (M) and WF' b (A) = 0, we have A e %°°' l (M) 

Note that WF b corresponds to the usual WF away from the boundary. 

Letting Lf(M) = [u G L 2 oc (X°); f x \u\ 2 ^-da], we have corresponding b- 
Sobolev spaces defined by 

Definition 3.3. u e H™' l (M) <=> V™- l (M)u C L 2 {M). 

We can define b-Wavefront set in b T*M for a distribution in an analogous 
way too: 

Definition 3.4. WF™' ! (n) = jp e b T*M\3A e ^°'°elliptic at p, Au e H™' l (M)} 
and the usual properties hold: 
(I) WF^' l (u) is closed and conic, 
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(II) WF^\u) n T*M° = WF^{u) 

(III) n m ,i K\M) = C~(M), U m , ; fl^M = C-~(M) 

(IV) A e vf< l '(M) maps A : # 6 m ' J (M) -> H™~ m ' (M) 

(V) for m < m' C WF^u) 

(VI) for u e H;°°' l \u) and A e V k > l (u) 

WF^' l+l '(Au) C n WFf +k ' l '(u), (34) 

and 

c (Ell b A) c . (35) 

(VII) if™'' If™ is compact if m < m' and / < /'. In particular, operators 
in W™' 1 are compact if m < 0, Z > 0. 

3.3 Domains 

We now establish a domain X> in which all of our arguments will take place. Let 
V denote the Friedrichs form domain of A + ^ffi , so the closure of C°° (M) with 
respect to the quadratic form: 

£(«):=/ |Vu(x)| 2 + ^^| M (x)| 2 da; (36) 

We then define X> by adding in D t : norm ||tt|||> = ||Aw|| 2 + ||Vu|| 2 2 + ||^m|| 2 2 - 
The s/2-th power of the operator A and D t will define V s analogously. 

Lemma 3.5. For u e C°°(M) and n > 3, we have 

^INI| < + ll-AHl 2 + \\D r u\\ 2 + ||Au|| 2 < C\\u\\l (37) 

This requires Hardy's Inequality which accounts for the n > 3 dimensional 
condition in our theorem. 

Lemma 3.6. Hardy's Inequality (from ^):For u G H l (R n ), n > 3, u(0) = 
we have: 




where the constant -k is the best possible. 
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Proof: See BH1 for the proof using integration by parts. □ 

We have as a special case from proposition 1 of HB-P-S-TZU that under our 
assumptions of dimension above on / in our potential which guarantee that for the 
Laplace-Beltrami operator on the sphere A\ + f(r, 6) + \{n) 2 (A(n) := is 
a positive operator on every sphere, we have: 

Proposition 3.7. Equivalence of Norms: HB-P-S-TZ\l There are constants c\ and 
c 2 such that: 

Ci||Vu(a;)|| 2 < Q(u) < c 2 \\Vu{x)\\ 2 (39) 

So the domain of our space under powers of P := A + is equivalent to 
the standard homogeneous Sobolev norm based on the powers of A. 

Proof. To show !3.7[ we can pick c 2 = ^1 + ^pjpr j since from Hardy's Inequality: 



Qiu) < J \Vu\ 2 + ^jf H 2 < ( 1 + ||Vu|& (40) 

To find ci (see statement of proposition), we integrate radially. Because of the 
positivity of A' + f(r, 6) + A 2 on every sphere, there is some 5 > such that: 



Q(u) = f°° f \d r u\ 2 + hv e u\ 2 + ^P-\u\ 2 dadr (41) 
Jo J\x\=r r 2 r 2 

> \V e u\ 2 + (f(r,9) + \ 2 )\u\ 2 dadr (42) 

JO r J\ x \=r 

r°° s 2 r 

> / — / \u\ 2 dadr (43) 



JO r2 J\x\=r 

5 2 \\-u\\ 2 L , (44) 
r 



Letting ci = &2 A ,,, we have using the above: 



Q{v)-c x \\Vufv = (1 - Cl )Q(u) + Cl [ l^l\ u \ 2 dx (45) 

r \v\ 2 

> / (- Cl /(r,0) + (l- Cl )5 2 )^d:r>O (46) 

JR n \ x \ 



□ 
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Proof. In |3.5[ clearly the only term that arises that is not included in the definition 
of the domain is ||-w|| which by Hardy's Inequality is bounded by ||ZVu|| and 
hence by the domain norm. □ 

Note that when extending this domain, we require an additional ||w||l2 term. 
Hardy's inequality will hold with this extra term if we decompose u by an ap- 
propriate cutoff function, supported at the origin with sufficiently small mass, i.e. 
u = ipu + (1 — ip)u. 

We also define a local version of the norm: 

\\u\\-,ioc = ||<HI- (47) 

for G (M) fixed to contain the region of interest. 

b-Pseudodifferential operators interact nicely with this domain, in particular, 
0-th order operators are still bounded and commute as expected with b-Differential 
operators which can be seen by using expanding the domain definition and using 
the usual theorems. This allows us to define b-Sobolev spaces with respect to this 
domain: 

Definition 3.8. For m > 0, let H^ b c be the subspace ofV of compactly supported 
u with Au G V for A elliptic over supp{u). Then loc is the subspace ofV of 
u such that for any <fi G C C (M), (fiu G HJ? . 

For m < 0, let are u G C~°°(M) of the form u = u\ + Au 2 with 

Ui,v>2 G T*i oc and A G \l/^ m (M). The norm is: 

II till #™ = mi {\\ui\\f, + \\u 2 \\f> u = ui + Au 2 } (48) 

T>,b,c 

with H% b loc the space ofu G C-°°{M) such that <pu G Hg b Jot <p G C C °°(M). 

Note that these definitions are independent of the A chosen. We also require 
wavefront set with respect to these Sobolev spaces: 

Definition 3.9. Let u G H% b [o Jor some s. For q E b T*M/ {0}, q £ WF^ (u) if 
there exists an A G \?j*(M) such that A is elliptic near q and Au G T>. 
Lemma 3.10. V = rH b (M). 

Proof. Use the definition of the norm of rH b (M) and the norm equivalence lemma. 

□ 

Corollary 3.11. Operators in are bounded on V and on H^for all m. 

Proof. Use the standard square root argument with respect to the Sobolev domain 
and note that this is equivalent to our new domain. □ 



17 



3.4 Second Formulation 



We can now reformulate the theorem on the blownup space in terms of b-wavefront 
set. 

The co- sphere bundle is a more appropriate place to state theorems about the 
propagation of WF(u) since they are already conic by definition and so propaga- 
tion in radial directions is hence only a matter of presence at all in a cone. To get 
an appropriate set of coordinates on the co-sphere bundle we pick r which is a 
degree 1 homogeneous function that does not vanish on our characteristic set, so 
we can renormalize fibers using it by letting | := -, £ := -, f := ±1 although 
only the first two variables will be used. 

The notion of b-elliptic regularity gives us that the solution to our equation is 
regular outside the characteristic set. This restriction corresponds to the restriction 
to unit speed geodesies from our original formulation. Note that □ := D*D t — 
(D* r D r - ±A 6 + where D* = D r + 2=± has a b-principal symbol of r 2 - 

£j — •1%-. Multiplying by r 2 and dividing by r gives the following characteristic 
set: 

S := {(r,^,t,e,0,r)| ( r 2 (r 2 n + /(r^))=0} (49) 
= {( r ,^,t,e,C,r)|r 2 -e 2 -|Clfc = 0} (50) 

where k denotes the length induced by the metric on the base § n_1 , so our so- 
lutions have WF b f ) {u) E S. On the interior of our blowup (away from the new 
boundary), we have WF b (u)n b T^ Q X = n(WF b f) (u)nT^ X) where vr is the blow- 
down map, so b-regularity translates into standard C°° regularity and vice-versa. 
The Duistermaat-Hormander theorem gives that WF b ^{u) is the usual maximally 
extended family of bicharacteristics inside £ on the interior of our blown up W 1 . 
It is useful to know an exact for for the Hamilton vector field here. 

Lemma 3.12. hM-V-WV If A G ^™(M) then the Hamilton vector field H a ofa = 
cr^ m (A) defined initially on the interior of the blown up cotangent bundle T*M° 
extends to the boundary as an element of Ub( b T* M) in coordinates (x, t, z, £, r, Q 
as: 

H a = (d(.a)rd r + {d T a)d t + (d (t a)d ei - (rd r a)d^ - {d t a)d T - (d 9j a)d (j (51) 
Proof Use the definition in coordinates and the b-projection map. □ 
So written as a b-operator, our Hamilton vector field is: 

H^a+i^l) ~ 2r(j t ~ ^ rd r - + ~2 °b ~2 °i ( 52 ) 
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Furthermore: 



Lemma 3.13. Angular momentum ( = for any bicharacteristics that strike the 
origin under this Hamilton vector field. 



Proof. To see this, lets analyze the flow for: 

2£ 



2rd t - -^rdr 



y Zi ry Zj 3 ry 

Call the propagating variable s and ' denote push-forward by d s . Our system is 



( t 1 


= T 


r' 


= _i 

r 




_ k^Q 


~ 2r 2 


t' 


= 




(5 2 +ICI 2 ) 


- r 2 





(54) 



We reparametrize flow with respect to the a parameter such that ds = r 2 dr, and 
rescaled this becomes: 

( t' = r 2 r 
= -r£ 



T' 

e 
l c 



o 

-e + \c\ 2 



(55) 



When > 0, the coupled system in r and £ produce: 

log(r(s)) = f£(s)ds 

m = >/ICUM>/ICU(*-ci)) 

which means 



Iog(r( a )) = tog(| 8 ec( ) /ici2( a - Cl ))|)-p 2 
Since sec is never 0,r ^ along this flow. 



(56) 

(57) 
□ 



There is a radial point at r = 0, ( = for where null-bicharacteristics strike 
the origin. This prevents us from using a diffeomorphism which orients the vector 
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field in some canonical direction and propagating through the origin as in the 
standard proof of Duistermaa-Hormander. 

We pick — £ = — ^ as our propagating variable, and note that H^ a+ f( r ,e)^(—£) 
is decreasing. This fact is crucial to constructing an appropriate commutant for 
our commutator argument, and it corresponds to the size of Ki $ ' o^ relative to 
the other terms. By picking the support of the commutant appropriately, we can 
show that this is the dominant term in our positive commutator argument. Also 
note that we're in the characteristic set along bicharacteristics with no angular 
momentum so r 2 = £ 2 . Notice that £ vanishes as we approach the origin instead 
of flipping and as we change sign, it accounts for outgoing bicharacteristics with 
the opposite momentum, and the absence of wavefront set where £ > will imply 
that there are no singularities on outgoing null-bicharacteristics. Restated again, 
our theorem becomes: 

Theorem 3.14. Let u in our domain be a solution of + j u = and 
over q = (t ,r ) let U be a neighborhood in S%(M. x [M n ,0]) around the set 
Q = { r = 0, t = t , r = ±1, £ = 0,0e S n ~\ C = 0}. Let U := U n E, £ the 
characteristic set defined before. In the region where t > 0: 

u n {r > 0, -| < 0} n WF(«) = =>• Q n WF^(m) = 0. (58) 

3.5 Proof of Second Diffractive Theorem 
3.5.1 overview 

We will inductively prove the reformulated theorem above inductively using the 
relative wavefront set. Because we are assuming all our solutions belong to V, 
they must be in some Sobolev space to begin with. Assuming (without loss of 
generality) that Q n WF° b {u) = 0, we want to show that Q n WF\ /2 {u) = 0. This 
would then guarantee inclusion in all orders of Sobolev regularity, showing the 
propagation of C°° regularity. 

Lemma 3.15. There exists A e ^°(M) such that: 

i([n + l±S-,A*A]u,u) = -(B*Bu,u) - ^(G*G>,u) (59) 

3 

+(E 1 u, u) + (E 2 u, u) + (Ru, u) (60) 
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In particular, \\Bu\\ detects order 1/2 wavefront set around Q and is bounded 
since every other term can be shown to be bounded. Therefore, we gain an addi- 
tional 1/2 order of regularity around Q. Note that the norm used here is relative 
to the domain mathcalD defined earlier. 

On the left hand side: assuming there is sufficient regularity in u for the left 
hand side to make sense, the left hand side will evaluate to by using the self 

adjointness of our operators and that (u + •^ 6> - > j u = 0. 

On the right hand side: we will construct this commutant using microlocal 
cutoffs, (i.e. one that cuts of in both spatial and phase variables). By its design, 
the commutant will expand into a number of terms on the right hand side. The first 

1 /2 

term —(B*Bu, u) is the main term which detects WF b (u) around the set Q; it is 
our goal to bound this main term and hence show an additional 1/2 order of regu- 
larity. This can be done by showing that all the other terms are either bounded or 
positive with the same sign as our main term. The second term — J2i(G*GiU, u) 
is a term with the same fixed sign as the main term so it helps our estimate. This 
term will be adjustable in size from how we construct A, and we absorb lower 
order terms from differentiation with non-propagating variables in the Hamilton 
vector field into this term. The third term (Eiu, u) will be supported on the in- 
coming co-tangent neighborhood and is bounded by our wavefront hypothesis 
U n {r > 0, -£ < 0} n WF(u) = 0. The fourth term u) is supported off 
of the characteristic set E and is bounded by elliptic regularity. The final term 
{Ru, u) is meant to be lower order, and ultimately bounded by inductive hypoth- 
esis. 

There will be technical caveats to consider regarding the final remainder term 
{Ru, u) which prevent us from working on the b-principal symbol level exclu- 
sively. However, suppressing these matters for the time being, we can provide a 
sketch of the lemma and construct A on the level of its b-principal symbol. 

First, let x( x ) be a smooth bump function with smooth square-root such that: 

*<*> = { J (-oo!'-2]U[2,oo) < 61 > 

where we call x'( x ) '■= 4>i{ x ) — 4>2 ( x )> eacri 4> smooth, and x( x ) be a smoothed 
Heaviside function 

*> = {o (-£'<.] < 62 > 

for which x'( x ) '■= 4>\{ x )i 03 a l so smooth. 
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In selecting x( x )> we require its derivative be the pair of bump functions with 
opposite signs x'( x ) = 4>\{ x ) ~ ^iK 3 -)' w ^ m 4>i( x ) supported on [—2,-1] and 
<f>2 (x) supported on [1,2]. Similarly, let x(x) and its derivative have a smooth 
square root supported on [0, 1] 

Taking the above remarks into account, we use the following construction for 
our commutant's b-principal symbol for some r > to ensure uniform bounded- 
ness of r from 0: 

a(r, z, t, £, C, r) := e^x (j j X {-r 2 + «£ + 25) (63) 
X (-(* - t ) 2 + a£ + 25) X (r - r )x ^ ~ ^~ \C\ 2 ^ ( 64 ) 

with A some constant. We abbreviate this by: 

a := e^x^XrXtXrXv (65) 

This symbol will localize in any given co-sphere neighborhood of our set of inter- 
est. It is clear that £ is localized by 5. The variables r and t — t are localized by 
the interaction of their parabolic cutoffs with £'s vertical cutoffs. The maximum 
width of the parabola cut off by the vertical is at most \/a5 + 25, which we can 
control for any given a. 

The first term, e cs - for C > a constant, decreases in — £ to ensure a main 
term with a consistent negative sign under differentiation by —d^ near our point 
of interest, and hence creates a nonzero symbol for detecting WF at our point 
of interest in our main term B. (Note that being in the characteristic set gives 
s + iClfc = r ' 2 so on bicharacteri sties that strike the origin, we localize in r as a 
result. ) 

The second term, x^ localizes in £, and under differentiation we get a term 
with support in the wavefront hypothesis region (left vertical edge of diagram) as 
well as one with the same sign as our main term (right vertical edge of diagram). 
The former is the symbol for E 1 and the latter contributes to the G. 

The terms, XrXt cut off in the slow variables, and the terms a£ introduce 
a parabolic tapering with respect to propagation in — £. This ensures that the 
Hamilton vector field points out of our support, and thus creates a bump function 
which is consistent in sign when hitting falling edge. These terms will make a 
small contribution to G, but will be dominated by the contribution from — £. 
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Figure 6: Top: tapering support ensures differentiation against a falling edge. 
(See end of chapter note for more detail about behavior in the r variable.) Mid- 
dle: profile of commutant. Bottom: profile of commutant under differentiation 
by Hamilton vector field; the bump on the left is supported where we have wave- 
front assumption, the non-zero middle portion detects wavefront set at the point 
of interest, the bump on the right has the same sign due as the middle portion. The 
like signed bump extends along the curved edges of the support (not pictured here) 
thanks to differentiation occurring only against the falling edge of the support. 
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The term x(t) keeps r > 0. Finally, xx, cuts off fiber variables to ensure that 
we have a symbol; the 5 parameter is meant to adjust the support of \ so that it 
is in some small tubular neighborhood of the characteristic set. Since this term 
is constant in a neighborhood near the characteristic set, differentiating it yields 
bumps supported off of the characteristic set whose quantizations are bounded by 
elliptic regularity. This term also keeps r > when — £ is negative and enables us 
to use the WF hypothesis on E 1 . 

We can directly see these terms appear by a direct application of H D+ / to the 
symbol we constructed: 



where: 



H D+ ^(a) = -b 2 - g + e x + e 2 



(66) 



b 2 := 
9 ■= 



ei := 
e 2 



-2 



e+\c\'- 



+Xi<f>l ( 




A 



T 



' 2 



'e+\c\'- 



r 2 + a£ + 25 



-2£ + 2 



XrXtXrXT, 

e + \cr~ 



a 



+XiXAl Ut-t ) 2 + aZ + 25 



a 



2(t-t )r + 2- 

T 




XrXtXrXY, 



e ^XtXrXtXrH D+ ^(xz) 



(67) 
(68) 

(69) 

XtXrX-B (70) 
(71) 



XtXt, 



(72) 

(73) 
(74) 



Note that for g, since r > 0, r — r2 < 5, and £ as well as t — t are cut 
off from being too big by Xr and xt respectively (in conjunction with xg), we can 
choose a to be as large as we'd like and force the entire term to be negative. 

The final term e 2 is clearly supported off of S since xs is constant in a tubular 
neighborhood of £ and hence any derivatives will only create terms off of E. 
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It is tempting to now quantize in the b-Calculus and finish the argument stated 
in lemma 13 . 1 51 however (Ru, u) requires additional care. This is because our 
domain involves b-differential operators of the form D t , -rD r , and -Dg^ but the 
lower order error term arising due to the quantization only a priori obey R e 
■%^ r ^ 1 (M) because of a deficiency in the filtration from the principal symbol map 
on the b-calculus. Using the operations of the calculus against a pair of operators 
will presume the "worse" weight of the two in the result, and this means that it is 
possible non domain bounded terms such as -D t Op(-) are present in R. 

There is a more technical argument on the sub-principal level we will now 
present which hinges on factorizing out the domain related D t , D r , and -Dg i 
terms at each step of the iterative argument. An elliptic regularity lemma from in 
HM-V-W1 lets us relate and these terms. In particular, we end up with a bound on 
a single term involving D t which is elliptic since |t| > and hence detects WF 
around the origin. 

The remainder also includes the inverse square potential adjoined to □, which 
can be absorbed into \rD r by an application of Hardy's Inequality as we shall see 
in the full proof. 

3.5.2 Setup For the Factorized of Main Propagation Theorem 13.141 

Our factorized argument will be similar to what was presented in the overview by 
iterating on differential order, but now we explicitly factor out D t , -rD r , and -Dg i 
at iterative step to ensure a domain bounded term for all lower order remainders. 
We denote these good differential operators by Q; L . 
In the end, we want to investigate the behavior of 

(-i[A*A,n + i^±]u,u) (75) 

iteratively while maintaining a factorization by Qi's. This involves expanding the 
operators on the right hand side by Leibnitz rule and investigating the properties 
of each resulting term. We further expand each term on the principal symbol level 
and classify these terms into the categories by the template: 

H n+ f (a) = -b 2 -g + ei + e 2 (76) 

seen as seen in the previous section. A lemma will allow us to conclude that 
the terms corresponding to G have a consistently "good sign" (i.e. are always 
positive). The terms E\ and E 2 are bounded by wavefront hypothesis and ellip- 
tic regularity respectively. Finally the R term is always domain bounded by the 
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way we iterated this argument only after factorizing out terms which are well be- 
haved in the domain. This shows that the elliptic B terms which will finally be 
{D t TBu, D t TBu) are bounded and hence detect a higher order of smoothness 
than before at the origin. 

We choose A e ^(M) with principal symbol a = cr ,o(^4) as before: 



a(r, z, t, £, C, t) := e^x (j j X (~r 2 + «| + 25) (77) 

X (-(* - to) 2 + a£ + 25) X (r - r ) X ^ ~ *]~ ^ ( 78 ) 

The b-quantization of this symbol requires us to use charts on S™" 1 . First we 
take {(f)i(x)} a partition of unity supported in these coordinate charts. Let {Xi( x )} 
be smooth cutoffs equal to 1 on a neighborhood of supp(0j) and still supported 
in the coordinate charts. Let p e C£°((— 1/2, 1/2)) identically 1 near 0. Now in 
local cooridnates, we can define the b-quantization using projective coordiantes 
in r: 

A, := (27r)-»y e i^<^h^p (^—^ a^U,^)^^)^ 

(79) 

with er some choice of right density. (Note that our a is independent of 9. We've 
also just blown up in spatial variables, not time.) Define globally 

Optrv(a) = ^XiAXi (80) 

i 

Recall that □ = D* t D t - D* r D r - 4, A§n-i where D* = D r + and 
Agn-i is the spherical Laplacian which in coordinates is of the form A s «-i = 
Y]- ■ —j=D e k lJ \/\k~\D e . for k tj G C°° on some chart and \k\ the determinant of 

k ij . 

We first investigate how our commutant commutes with any of our Q l, s. For 
A := Optrv{a) as before, we now consider A* A (which is very similar on the 
principal symbol level to A) for self-adjointness considerations, and by Leibnitz 
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rule: 




(81) 



i[A*A, D* t ]D t + D*i[A*A, D t ] 
-i[A*A, D*]D r - D*i[A*A, D r ] 



(82) 
(83) 






k\D .]+i[A*AJ(r,9)} 



k\D ej +f(r,9) 



) 



(85) 



(84) 



We will need some lemmas about the behavior of the [A* A, Q % \ terms present 
above. 

The following lemma about the commutators involving D r from HM-V-WB 
allow us to improve on the order of the weight more than just what the b-calculus' 
filtration tells us. Note that as a b-operator D r = \rD r and so a priori in the 
b-Calculus filtration [D r , A] e ^"^(M). 



Lemma 3.16. hM-V-WJ IfA e ^(M), there exist B e ^(M), C E ^~ X {M) 



It is crucial that B and C on the right hand side are weightless operators. The 
proof of this lemma makes use of the fact that [rD r , A] has an additional order of 
weight regularity than the b-Calculus would a priori give us. This lemma is shown 
by using normal operators which have their coefficients involving the r variable 
are frozen at r = 0. Please refer to HM-V-WI for further details. 

The next lemma deals with the behavior of [A, Dg}. Our commutant A has 
principal symbol that is constant in the spherical variables Of, such symbols are 
called basic symbols so this means that Dq.(o) = on fiber variables on each 
chart. We are still left with some non characteristic terms though: 



depending continuously on A such that 



[A, D r ] =B + CD r 



(86) 



with a(B) = -id r (a(A)), a{C) = -id^(a(A)). 



Proof. See HM-V-WD for the proof involving normal operators. 



□ 
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Lemma 3.17. HM-V-WV IfA G ^Sf(M) where A = Op(aip ( r2 ~^ +l ^ 2 )) where a 
is independent of 9 and tp cuts off near S and 

[A,D e ^] = Be^T(M) (87) 

oZ>eys WF(S) n E = 

Proof. Recall the quantization on each chart: 

(88) 

and 

Optrv(a) = ^XiAXi (89) 

i 

Note the independence of a from 9 and £, and so such operators commute with 
Dg at the cost of terms [%, A?] which are smoothing terms because the Schwartz 
Kernel is smooth. These will not affect our lemma as they are residual. 

On each coordinate patch near the diagonal, we are able to represent the quan- 
tization under a change of variable in 9 as: 

A := (27r)-"| e <^-(-'))+^') p fr-±\ ~a(r,t,OMr,e,Z,(ndCda 

(90) 

Note our a has the form a : = a%s and that ipi = (fiiXz here has a change of 
variables under z. The sums of the derivatives on if) will cancel under summation. 
The only possible contributions which come from [A, Dg] must be from d{ip) 
where xt. is non-zero, which must be non-characteristic. □ 

A final lemma gives the behavior of [A, ^] 

Lemma 3.18. If A G ^(M), then: 

[A*A,±] = C (91) 
where C G ^ m ~\ a(C) = -^2id^(a(A)). 

Proof. We use the same lemma as in the proof of (|3.16l) . which specifies that 
[rD r , A] has an additional order of weight regularity than the b-Calculus would a 
priori give us. Here: 

[A,r 2 i] = [Ar 2 ]l + r 2 [A^] (92) 
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and so 

[A^] = ^4,r 2 ]^ (93) 

where the middle term on the right gains us 2 degree of regularity as before. The 
principal symbol comes immediately from lemma I3T21 □ 

Lemma 3.19. For the specific A* A G which we mentioned before, this yields: 

Principal Symbol Order 

[A* A, D*\ & [A* A, D t \ 

= U a{Li) = 2iad t {a) G *g(Af) 

[A* A, D*] & [A* A, D r ] 

= Bi + CiDP a(Bi) = 2iad r (a) B { G ^g(M) 

o-(C0 = 2iadt(a) d e ^(M) 

[A* A, 

= £Os <t(C 5 ) = 4za^(a) G ^(M) 

= E 6 + R WF(£ 6 ) n E = ffom'c operator) 

a(R) G *^(M) 

[i4M,/(r,fl)] 

= C 7 + E 7 a{d) = 2iad^a)rd r (f) G ^ b \M) 

WF(E 7 )nS = G^(M) 

(94) 



Proof. See the preceding lemmas for the first three rows. 

The non-characteristic terms are from d( terms which are supported outside 
our characteristic set as our symbol is basic. Note also that depend only 
on 6i's and the principal symbol of [A* A, k ij ], [A* A, y/\k\], and [A* A, -±=] axe 



where Op^ = -j=D* 9i W ^J\k\D e 



2iad^(a)de(k lJ ), 2iad((a)dg(^\k\), and 

2iad^(a)de(l/ y/\k\) respectively. Since is a metric, its inverse is smooth and 
each of terms with derivatives in 9 will be a smooth function in 9 (with respect to 
trivialization on chart coordinates which we have built into our quantization). Our 
basic symbol then forces the highest order of these terms off the characteristic set. 
The final row comes immediately from 13 .121 □ 

Now we can write expressions for these factorized forms whose remainders 
only pick up additional domain bounded terms. Temporarily renaming C 3 and C 5 
to C 3 and C 5 reserve them for after we commute, we get: 
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-i[A*A,U+ f -^P-}= (95) 
[L x D t + D* t L 2 - {(B 3 + C 3 D* r )D r + D* r (B 4 + C 4 D* r ) + (96) 
^0>E £>;^D Zj - /(r, ^)) + E 6 + R}} + 1(C 7 + £ 7 ) (97) 

We can regroup to fit our final commutator argument better by letting E 3 = 
[D*, C 3 ] G ^ X (M) and £ 5 = [D* z ., C 5 ] G where the latter operator is 

supported off the characteristic set and so is absorbable in E e ; call C 3 = C s + i£ 3 
which we can also reabsorb by using that C is basic. This gives: 

-j[i*i,D + ^f] = [LiD t + DIL 2 — B 3 D r — D* r B± (98) 



-{D* r (C 3 + C A )D r + ~D* ei C 5 k ij ~D 6j (99) 



r 



f(r 0)1 1 
-^^C 5 -) + -(C 6 - E 6 + E 7 ) + fl](100) 

Here i? is a sum of terms from Q 1 ^ 1 and (Q*) 2 \I/~ 2 which is lower order and 
comes from the commutators of our explicitly factorized terms terms above. (In 
particular, we note that these are bounded with respect to lower order domain 
norm when paired as (Ru, u).) When we ultimately pair with our solution u, we 
want to absorb the terms with only first order differential operators in line (1981) 
into line (|99l) . We will see that this is because the principal symbol of the second 
line will dominate the smaller symbol of the first as in the principal level overview 
argument earlier. 

Let's examine the terms dominant C 3 + C 4 and C 5 present in (|98l) and see what 
contributions their (common) principal symbol Aiad^a makes. 
Recall our commutant: 

a(r, z, t, e, C, r) := e c * X ( f \ x[-r 2 + a£ + 25) (101) 



5 



x(~(t- to? + ai + 26) x(r - r ) X ( - ^ ^ ) (102) 
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abbreviated as: a := e c ^x^XrXtXrXT:- First we make a few factorizations to 
isolate the types of terms we are interested in: 



A 



Aiad^a = 4:t-e 2M xjx 2 r XtX 2 T xl 



T 



8t 



~2~2~2 2 
'X.r'X.t XtXy^ 



+4ie ci x}Xr(p3, 



~2 ~2 2 
XtXrXx 



~2 2 

X T Xz 



(103) 

(104) 

(105) 
(106) 
(107) 



+&e c t X \x 2 rXtfo,t 

+Aie c ^x\x 2 rX 2 tXT.X% (Xs) 

Recalling the template for the terms we want in our pairing from the lemma in the 
previous section, this suggests we decompose each into the following: 

_ iC. = —T^B* BT-i -G + E 1 + E 2 + R (108) 

The term T_i G ^\M) with cr(T_i) = ^Xs 

(the xy, makes this a symbol) 
is elliptic. We will use this term in subsequent sections to shift orders on paired 
terms. 

1 12 

As before, B E $ b (M) is the main term quantized from 



b := VrAe M x^XrXtXrXs 



(109) 



i.e. a(B*B) is — it 2 times line (11031) and is elliptic at the boundary r = 0(from 
our assumption that r > 0) so that it detects 1/2 order of extra regularity there. 

G E $~ l (M) the "good sign term" from (11041 ) (11051 ) and (I106I ). It corresponds 
only to the terms arising from the £ derivative, and we define its principal symbol 
#by: 



9 :z 



-Mi 



Xz<Pl 



5r 



~2~2~2 2 
XrXt XtX^ £ 



x^xAi (■) 

2A£ 2~2~ j 2 



a 

it j 



~2 ~2 2 

Xt XtX^q 



a 
it 



~2 2 
XtXc^ 



Ei E % X {M) the ll WF hypothesis term" from (11041) . 



2 

57 



e 2 + icp 



~2 ~2 ~2 2 



(110) 

(111) 

(112) 
(113) 
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E 2 G ty b 1 (M) from (11071) the non characteristic term: 



C£ 2 ~2 ~2 ~2 q / \ 



(114) 



andi? e \T/^ 2 (M) a domain bounded lower order remainder. (See previous section 
for a diagram.) 

The main summand above involving C3 + C4 and C5 can be rearranged at the 
cost of commutator terms to form: 



where Kij correspond to the smooth functions that appear in the Laplacian as be- 
fore. We group together terms with properties corresponding to E 1 and E 2 . We 
have also grouped resulting commutator terms (which are all of small enough 
order) into the remainder term, which are now are sums of the remainder terms 
before and also new terms R £ Q*Qi^l(M). (These are also bounded by our 
domain norm) Note that in particular, B is basic (its symbol constant on fiber 
variables) so commutators with Dg i are off the characteristic set and we can com- 
bine these into our existing non characteristic set term. 

Note that all subsequent terms named E\, E 2 ,R will not necessarily the same 
expression in every line, they will simply obey the corresponding properties cor- 
responding to terms bounded by wavefront hypothesis, bounded by elliptic reg- 
ularity (from having support off £), and being a lower order domain bounded 
remainder term respectively. 

Since we want our final bound to be on the term (D t TBu, D t TBu), we want 
to exchange the differential operators D r and \D Qi in line (1981) to D t . This can be 



done using ( □ + 2 - ) u = to help us. As we know the explicit form of the 




(115) 



(116) 




Laplacian - (□ + ^) = D* t D t - D* r D r - \D e W\D d] + 



/M) 



r 



, we can 



exchange D r and ^Dg terms for D t by writing the above as: 




(117) 



G (D*D t ) + E l + E 2 + R 



(118) 
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where W is given by B*T* l T_iB + G where we have used the fact that B is basic 
and T ! is lower order to absorb commutators from commuting terms around into 
the E and R terms again. 

Likewise, Li and Bi have symbols: 

2iad t a = 2ie 2M X \xlMi,t{-)W--to]xlx 2 a (119) 
2iad r a = 2te 2Ai X 2 i XrM)[^}XtXlx 2 a (120) 

Note the similarity to the terms which we termed G t and G r above, and in final 
pairing expressions, these terms will be absorbed into the G terms which have 
good sign and make no contribution against our main commutator. 

3.5.3 Main argument for the Propagation Theorem 13.141 

Now for the main inductive argument where we introduce the pairing with u, a 
domain bounded solution. For now assume s = without loss of generality. (We 
will shift s with a separate operator at the end which smooths as well.) Because 
of the restriction of our solution to a domain, we can assume that for some small 
s , that at our point q WF S ^. Now we assume WF s b ^u PI U = 0. We will 

show that WFf~£ 2 w D U = for some open neighborhood U C U Let's pair the 
actual commutator expression we will use with our solution, expanding using the 
expressions we established above: 

-i([AM 1 D + M] lt)U ) = (121) 

(-W(n + i^l)u,u) (122) 

— {D t T-xBu, D t T^Bu) (123) 
-{GD*D t u,u) (124) 
+(LiD t u, u) + (L 2 D t u, u) + (B 3 D r u, u) + (B 4 u, D r u) (125) 

+ (^C 7 u,u) (126) 

+{Etu, u) + (E 2 u, u) + (-jEfu, u) + (Ru, u) {111) 

Note that this is only a formal expression as we don't necessarily have suffi- 
ciently regularity on u to know that these L? pairings make sense; we will address 
this issue with an additional regularizing term and a weak convergence argument 
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at the end of the proof. Here we have absorbed the non-characteristic E$ and E 7 
term into E 2 . (The weight of \ will make no difference as after using the cauchy- 
schwarz followed by Hardy's inequality, we can substitute - with D r plus some 
multiple of the domain norm of u which is inductively bounded.) 

That G be positive is crucial. Using the form of the principal symbol from 
( II 101) , we decompose G into 3 operators with positive principal symbol G = 
G^G^ + G* r G r + G* t G t , where G^, G r , and G t and their adjoints are quantized 
from the square roots of —i times the principal symbols corresponding to the lines 
starting at (II 101) . We will ultimately absorb the LiD t and BiD r terms from line 
(11251) into the G t and G r terms respectively. By showing LiD t and BiD r terms 
are smaller in norm than the norm of the corresponding positive G t D t and G r D t 
terms respectively in (11241) . we see that their sum with G is still positive. 

The naive, purely principal symbol level argument works for the D t terms as 
it doesn't have a b-weight. First consider the terms (LiD t u, u) which we now 
absorb into G t . To do this, we consider the difference \D* t G* t G t D t — LiD t as 
there are two terms and perform one step of a square root argument to achieve a 
positive term which helps us as it will have the same sign as G t plus a residual 
term with lower order that has no weight and is domain bounded by inductive 
hypothesis. We get 

((D* t G* t G t - Li)D t u, u) = (S*Su, u) + (Ru, u) (128) 

where *{S) = ^a{{G t D$D t - L t D t )) = y^xM (■) [2a - 2(t - t )]x?x|^, 

R G D t ^ 1 (M). This is a smooth symbol as r > and we can make a as large 
as we'd like to force the factor 2a — 2(t — to) to be positive. Therefore, this term 
will have the same sign as our main term and the residual term is bounded by 
hypothesis. The same technique works for the other term involving D t . 

We could try to use the same technique for the terms (BiD r u,u), however 
as a b- operator, D r = \rD r and we would impose an arbitrary weight on our 
remainder which is not attached to the specific differential operators mentioned 
above, which was the whole point of factorizing Q % at each step. Instead, we 
show an estimate based on a lemma that allows us to estimate domain derivative 
in terms of each other, in particular D r and -Dg i can be bounded interms of D t . 

We factorize Bi in terms of G which becomes: 

B, = G*T^RiT x G r + R (129) 

where cr(Ri) = ^Xs, 0"(T X ) = (r)xs and R G ^^(M) is residual. Clearly the 
supremum of the symbol of cr(Ri) = ^Xs i s arbitrarily small since r is controlled 
by Xr' s support, and we can make a as large as we'd like. 
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Letting v = G r T\ and commuting D r using lemma I3TT61 and absorbing these 
weightless lower order terms into R yields: 



(BiD r u,u) = (RiD r v,v) + (Ru,u) (130) 

< Ili^A^HIMI (131) 

by cauchy-schwarz. 

Ti is elliptic so it has a parametrix T"_i and Id = T{T_i + F. Using this and 
letting v = G r Tiu: 

\\RiD r v\\ (132) 

= \\R i (T 1 T. 1 + F)D r v\\ (133) 

< \\(RiTi)(T^iD r v\\ + \\RiFD r v\\ (134) 

< 2sup|o-(i2 i )||||T_iD r v|| + ||i^r_iD r t;|| + ||i2 i i5 , D r t;|| (135) 

for somei?^ G %\M). 

A useful proposition for cross-terms \ \v\\ is: 

Proposition 3.20. ||w||||f|| < ^IM| 2 + 7lM| 2 

Proof. Expand: (-^Hl - v^lMl) 2 > 0. □ 
Applying this to \\RiD r v \\ \\v\\: 



(RiD r v,v) < 2s\xp\cx(Ri)\\T- 1 D r v\\\\v\\ + 2 7 |H| 2 (136) 

+-\\R' i T- 1 D r v\\ 2 + -\\FD r v\\ 2 (137) 

7 7 

We still want to absorb this term into ||Gf-Diw|| 2 , and this will require that we 
replace D r in our expression by D t . This requires the following crucial lemma 
which lets us do this without introducing non-domain bounded remainders: 

Lemma 3.21. M-V-W\l For K C U C b S*M, K compact and U open and 
A r G ^> s b ~ l (M) for r G (0, 1] a basic family with WF' b (A r ) C k bounded in ^ s boo . 

Then there exists Gg$j 1 ^ 2 (M) with WFj,(G) C U and Cq > such that when 
Uu = 0, 

/ (\d M A r u\ 2 - \D t A r u\ 2 ) 

J M 



<C Q {\\u\\l ioc +\\Gu\\l) (138) 
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Proof. For a proof, please refer to HM-V-WL Note that the r parameter is required 
for a later regularity approximation argument. □ 



This allows us to estimate the D r term as follows: 



/ (\d M A r u\ 2 - \\D t A r u\\ 2 + \\D t A r u\ 
Jm 

|2 , ||i^f..||2 \ , II n A ..112 



(139) 



< Co(lkll4 oc + H Gu lfe) + II D t A ru\r (140) 

From the expression (11361) . only the first two terms aren't bounded by hypothesis 
on their order: 

2sup\cr{R i )\\T_ 1 D r G r T 1 u\\\\G r T 1 u\\ + 2j\\G r T 1 u\\ 2 (141) 

In the first term we can commute D r using our lemma I3T61 at the cost of a resid- 
ual term which can be bounded by background regularity. We note the that first 
term is a cross term again so we may use 13.201 as before to express this in terms 
of squared norms. Then after using the lemma, the expression is bounded by 
2Csup |cr(-Rj)|||G>.D t w|| 2 plus lower order terms. We are left with a 7||G r T 1 M|| 2 
term which group with the second term. 

In the second term, we can commute as desired and we note that T\ has a 
symbol that asymptotically matches D t . Thus, 27||G>T 1 ii|| 2 can be expressed as 
27 1| G>.D t M || 2 with a residual term that is bounded. 

The term 2C sup \a(Ri)\ \\G r D t u\\ 2 can be made arbitrarily small by 2 sup \ a(Ri 
which becomes vanishingly small for large a and by 7 which can clearly be 
as small as we want in 13.201 They can therefore be absorbed into the existing 
|| G>.D t w || 2 term with coefficient 1 without changing its sign. 

The other term involving B A and D r is handled similarly. 

This leaves the term {^Cju.u) from our potential. This may seem lower 
order, but it is not sufficiently low to be absorbed into our domain norm under the 
wavefront assumption we had. (That would require C 7 G \1/^ 2 (M).) To absorb 
this term into existing terms, we factor C 7 into the same terms as we did for 
with the symbol Aiad^(a): 

C 7 = —RT* 1 B*BT_i - RG*G + E l + E 2 + R' (142) 

with a(R) = rd r f(r,6) All terms except RT*B*BT and RG*G end up being 
bounded for the assumptions stated before. G decomposes into three terms as 
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before G*Gi- After commuting (without adding weights), by cauchy-schwarz 
we can bound the pairings by: 

(-RBT-xUi-BT^tu) < W-RBT-xuWW-BT-xuW (143) 

ry ry ' ry ry 

(-RGu,-Gu) < \\-RGu\\\\-Gu\\ (144) 

ry ry ry ry 

and we would like show these terms are arbitrarily small in order to absorb them 
into the main terms in the pairing. We can bound these terms in the pairing in 
terms of squared norms using 13.201 as before. Then, by Hardy's inequality, each 
of these norm expressions are bounded by 

7 2 (\\D r BT_ lU \\ 2 + \\BT^u\\ 2 ) (145) 
7 2 (\D r G r u\\ 2 + \\BT-mW 2 ) (146) 



We may exchange the first terms in each expression above for terms involving 
D t and a lower order term finite in the domain using the lemma above, so they 
become 

7 2 (|| D t BT_ lU \\ 2 + \\BT-xuW 2 ) + Lo -t- ( 147 ) 
7 2 (j| D t Gu\\ 2 + \\BT_ lU \\ 2 ^ + l.o.t. (148) 

respectively. Making 7 small controls these terms. 

After applying |3 .2 1[ ours is a question of making the following terms arbitrar- 
ily small for potentially large \ values: 



7 
1 

7 



2U D t RBT_ lU \\ 2 (149) 

\\\D t RGu\\ 2 (150) 

We would like to absorb these into corresponding terms involving BD t and GD t 
in our main expression whose signs we must maintain in order to preserve signs. 
The rest of the terms are bounded by being weightless and lower in order. 

These will indeed be small since the supremum of cr(R) can be made arbitrar- 
ily small by our control in r. With out loss of generality, we consider the term 
involving G after (weightless) commuting as usual. 
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Lemma 3.22. For R G we have M s G *°(M) -R', i?" G 1 f/ie inequality: 
\\RD t Gv\\<2 sup |CT(i?)|||AG^|| + H^AftfAGull + H^'A^II (151) 

supp(cr(B)) 

Proof. To see this, create Ms = Op(xr"') G \I/°(M) using smooth cutoffs 
that equal 1 on the support of cr(G) so that a(M)a(G) = a(G). (We could for 
instance double the width of the domain profile for each x-) ana i n particular we 
are concerned with cutting off in r to bound the size of rd r (a). Then MG = 
G — R' for some R' = ^ 3 ' 2 (M). Thus we can express the above as: 

\\RD t Gv\\ (152) 

< \\RD t MGv\\ + \\RD t R'v\\ (153) 

< sup \a(R)\\\\D t Gv\\ + \\R"D t MGv\\ + \\RD t R'v\\ (154) 

supp{a(G)) 

for some R" G \I/^ 1 (M). □ 

The coefficient sup supp{a{d)) cr(R) = ^(swp r£mpp{a{M)) 5d r f(r,9)) is clearly 
arbitrarily small by making 5 small as r is localized by the principal symbol of G. 
The dangling terms ||_RGu|| \\RD t R'u\\ and \\R" D t MBT-iu\\ are of lower order 
and hence bounded. 

Squaring \\RD t BT_iv\\ and HiiZ^GvH as presented in our lemma yields our 
main term and "good sign term" with arbitrarily small coefficients as well as cross 
terms with of the aforementioned with bounded terms. The cross terms can again 
be bounded by 13.201 Since the resulting main and "good sign" term can be made 
arbitrarily small, we can absorb them into the existing main and "good sign" terms 
in the main argument while maintaining their signs. (All absorption steps are 
clearly compatible with all other absorptions since we could pick our symbols 
arbitrarily small for each.) 

Returning to our main pairing: 

-i([An+^K«) = (-W(a + £^)u,u) (155) 

— {D t TBu, D t TBu) (156) 
-Cj2{G*i D tu,GiD t u) (157) 
+ (E 1 u, u) + (E 2 u, u) + (Ru, u) (158) 

where CJ2(G*D t u, GiD t u) is positive from our argument just now. The expres- 
sion: (D t TBu, D t TBu) is thus bounded by the absolute values of the remaining 
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terms which are finite. A final application of the lemma above allows us to bound 
||d]RnTSu|| by \\D t TBu\\, and therefore the entire domain norm with 1/2 higher 
order. However, the pairings may not always make sense for insufficiently regular 
u. 

We finally address the required regularity of u for these pairings to make sense. 
We resolve this by using an approximating argument using sequence of smoothing 
operators A r quantized from the symbol: 

|r| s (l + a|r| 2 )- s / 2 (159) 

Letting A r = AA r : 

-i([A* r A r , □ + ^P-\u, u) = (160) 

(-W(D + ^^)A r u,A r u) (161) 

- (D t TBA r u, D t TBA r u) ( 1 62) 

-J2(G*D t A r u,GD t A r u} (163) 

+ (E 1 A r u, A r u) + (E 2 A r u, A r .u) + (RA r u, A r u) (164) 
The left hand side is: 

i([A* r A r ,n+^-^]u,u) (165) 
= (A*,A r (d + ^P-^J u,u)-{(d+ ^P^j A*A r u, u) (166) 

= (A r (n + ^^^ u ,A rU )-{A rU ,A r (n+ii^ju) (167) 

which is as u is a solution. 

The "good sign" term has the same sign as our main term so it does not hurt 
our inequality and we can ignore it. The other terms are bounded because of WF 
hypothesis in the case of E\, elliptic regularity in the case of E 2 , and inductive 
hypothesis for R. Therefore, the L 2 norm \\D t TBA r u\\ is uniformly bounded as 
r — > 0. This means that weak compactness holds and that some subsequence of 
D t TBA r u converge to weakly in L 2 , or strongly to some distribution. However, 
the Riesz representation theorem guarantees that this is some L 2 function, and by 
uniqueness of distributions, we have that D t TBA u is an L 2 function. 

Hence || D t TBu\\ 2 is bounded and we have no wavefront set of order \ +s. We 
can inductively repeat this argument by shrinking the elliptic set of B arbitrarily 
much at each step and conclude q WF^^(u). □ 
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